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We present results from a study of meson spectra and 
structure in the limit where one quark is infinitely heavy. The 
calculations, based on the framework of light-front QCD for- 
mulated on a transverse lattice, are the first non-perturbative 
studies of B-mesons in light-front QCD. We calculate the 
Isgur-Wise form factor, light-cone distribution amplitude, the 
heavy-quark parton distribution function and the leptonic de- 
cay constant of B-mesons. 



I. INTRODUCTION 

Heavy flavor physics has become one of our most 
promising tools in searching for physics beyond the stan- 
dard model. However, despite the use of heavy quark 
symmetry and heavy quark effective field theory, many 
uncertainties in extracting standard model parameters 
from data on heavy meson decays arise from the often 
unknown nonperturbative hadronic matrix elements for 
these decays. For example, in order to extract the CKM 
matrix matrix element \Vbc\ from semileptonic B D* 
transitions, it is necessary to extrapolate to the point 
v-v' = 1. Even in the heavy quark limit, this involves the 
Isgur-Wise form factor, which is not determined by sym- 
metry. For this and similar applications, it is essential 
to have reliable predictions for hadronic matrix elements 
that multiply interesting standard model parameters in 
the physical amplitudes.^ 

Euclidean lattice gauge theories have played an im- 
portant role for calculating some of the relevant non- 
perturbative matrix elements. However, especially for 
matrix elements that involve a large velocity transfer 
as well as matrix elements that probe light-like correla- 
tion functions (such as the B-meson distribution ampli- 
tude or the parton distribution functions, which are rele- 
vant for inclusive semileptonic decays), it would be much 
more natural to describe the nonperturbative hadronic 
physics in the light-front (LF) framework (or infinite mo- 
mentum frame), i.e. in a Hamiltonian framework where 
= + x'^ is 'time'. 

What distinguishes the LF framework from all other 
formulation of QCD is that the observables, which are 
dominated by light-like correlation functions, become 
equal 'time' correlation functions and hence have a very 



direct and physical connection to the microscopic degrees 
of freedom in terms of which the Hamiltonian is con- 
structed Q. Because of this unique feature, it should 
be much easier in this framework to gain a physical un- 
derstanding between experiment and phenomenology on 
the one hand and the underlying QCD dynamics on the 
other. 

Of course, just like in any other approach to QCD, it 
is necessary to regularize both UV and IR divergences 
before one can even attempt to perform nonperturba- 
tive calculations. The transverse lattice [|j is an attempt 
to combine advantages of the LF and lattice formula- 
tions of QCD. In this approach to LF-QCD the time and 
one space direction (say x^) are kept continuous, while 
the two 'transverse' directions x_l = {x^,x'^) are dis- 
cretized. Keeping the time and x^ directions continuous 
has the advantage of preserving manifest boost invari- 
ance for boosts in the x^ direction. Furthermore, since 
x^ = x^ ± x^ also remain continuous, this formulation 
still allows a canonical LF Hamiltonian approach. On 
the other hand, working on a position space lattice in 
the transverse direction allows one to introduce a gauge 
invariant cutoff on _L momenta — in a manner that is 
similar to Euclidean or Hamiltonian lattice gauge the- 
ory. 

In summary, the LF formulation has the advantage of 
utilizing degrees of freedom that are very physical since 
many high-energy scattering observables (such as deep- 
inelastic scattering cross sections) have very simple and 
intuitive interpretations as equal LF-time {x'^) correla- 
tion functions. Using a gauge invariant (position space-) 
lattice cutoff in the _L direction within the LF framework 
has the advantage of being able to avoid the notorious 
l/k~^ divergences from the gauge field in LF-gauge which 
lague many other Hamiltonian LF approaches to QCD 



^For a recent review, see Ref. 



The hybrid (continuous along with discrete) treatment 
of the longitudinal/transverse directions implies an anal- 
ogous hybrid treatment of the longitudinal versus trans- 
verse gauge field: the longitudinal gauge field degrees 
of freedom are the non-compact while the transverse 
gauge degrees of freedom are compact link-fields. Each of 
these degrees of freedom depend on two continuous {x^ ) 
and two discrete (n^) space-time variables, i.e. from a 
formal point of view the canonical transverse lattice for- 
mulation is equivalent to a large number of coupled 1 + 1 
dimensional gauge theories (the longitudinal gauge fields 
at each nj^) coupled to nonlinear a model degrees of free- 
dom (the link fields) |^. 
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This formaUsm has been previously applied to pure- 
glue QCD (numerical results for the static quark- 
antiquark potential can for example be found in Ref. |^] 
and for glueball spectra in Ref. |^). More recent applica- 
tions of the color dielectric transverse lattice formulation 
to mesons at large Nc can be found in Refs. where 
more details regarding the transverse lattice approach ap- 
proach can be found. 

In this paper, we make the first attempt to apply the 
_L lattice, used in the above light meson calculations, 
to mesons where one of the quarks is infinitely heavy. 
The paper is organized as follows. We briefly describe 
the Hamiltonian, the approximations used to solve it in 
section 2. Numerical results for a number of important 
hadronic matrix elements are presented in section 3. 



II. THE HAMILTONIAN 

The transverse lattice Hamiltonian for Wilson fermions 
that forms the basis of our work was first constructed in 
Ref. [ p2[ and was previously applied to mesons containing 
only light quarks in Refs. and we refer the reader to 

these works for more details. Here we restrict ourselves to 
a more qualitative description of the degrees of freedom 
that enter the Hamiltonian and some of the details of 
how heavy quarks are incorporated in this framework. 

The degrees of freedom that enter the Hamiltonian are 
quark and antiquark creation operators at each site and 
link-fields (which, in the spirit of the color-dielectric ap- 
proach we take as general matrix fields) on the trans- 
verse links connecting the sites. For simplicity, we work 
in the large Nc limit. The first approximation that we 
use is the same light front Tamm-Dancoff truncation of 
the Fock state expansion that was used in Refs. |p|,pl|, 
namely, a truncation upto and including the 3-particle 
sector. After reminding oneself that the fermion degrees 
of freedom (quarks and antiquarks) always occupy the 
lattice (transverse) sites and the gauge degrees of free- 
dom (link fields) reside on the lattice spacings, it is easy 
to understand the following two direct consequences of 
the aforementioned truncation: 

• The 2-particle states consist of a quark and an an- 
tiquark sitting on the same lattice site. 

• The 3-particle states consist of a quark and an an- 
tiquark separated by at most one link field. 

The Hamiltonian for the light quarks and link fields in 
this model contains the following 6 parameters 

• kinetic masses for the light quarks in the 2 and 3 
particle Fock sector 

• helicity flip and non-fiip hopping terms for the light 
quarks 



• kinetic mass for the link field 

which were determined in Ref. by imposing rotational 
invariance, fixing the physical string tension in _L lattice 
units in both the longitudinal and transverse direction 
and the tt and p masses. 

For the heavy quarks one has a variety of options. Con- 
ceptually, it would be cleanest to treat them as infinitely 
heavy fixed sources, using the fixed source formalism out- 
lined in Ref. 0. In such an approach, no new parameters 
would enter the Hamiltonian. In this work, we adopted 
a slightly different procedure to implement the heavy 
quarks: we treated them as finite mass quarks, but did 
not allow the heavy quarks to propagate in the trans- 
verse direction — only propagation in the longitudinal 
direction is allowed. The advantage if this hybrid proce- 
dure is that we could utilize computer code that was used 
to describe light-light mesons and the heavy quark limit 
results were then obtained by extrapolation. ^ Treat- 
ing heavy quarks consistently as finite quark mass heavy 
quarks would have the advantage that one could study 
1/Aff, corrections to the heavy quark limit. However, 
since this would not only require developing an analog of 
the NRQCD procedure on the LF, but also a determina- 
tion of a whole new set of parameters (the coefficients of 
the hopping terms for the heavy quark in the Hamilto- 
nian), we preferred to focus only on the oo limit, 
where no new parameters are introduced. 



III. NUMERICAL RESULTS 

The observables that we studied are the Isgur Wise 
form factor, the decay constant /s, the B-meson light 
cone distribution amplitude and the (light-quark) parton 
distribution function in a _B-meson. For all these observ- 
ables, we performed calculations at finite kinetic mass Mb 
(but no _L hopping) and then extrapolated to Mb — > oo. 
Any error bars/bands that we show for these observables 
reflect uncertainties arising from this extrapolation. 

In each of these calculations, we diagonalized the _L 
lattice Hamiltonian in a truncated basis, using continu- 
ous basis functions. As output of these calculations we 
obtained the wave functions (Fock space amplitudes) in 
the two and three particle Fock components tps{x) and 
ij)s{x,y) where the arguments in the two particle Fock 
component refer to the momentum fraction carried by 
the light quark and in the three particle Fock component 
refers to the momentum fractions carried by the light and 



• longitudinal gauge coupling 



■^Note that although this procedure produces the correct in- 
finite quark mass results through extrapolation, we cannot 
use the results to study finite Mb corrections, since the dy- 
namics of finite Mt, is not properly described (no transverse 
hopping!). 
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heavy quark. The discrete index s labels the spin con- 
figurations (4 in the two particle Fock component) and 
the spin/orientation components (4 x 4 = 16 in the three 
particle Fock component). 

In the following, we will use numerical results for these 
wave functions to calculate some important hadronic ma- 
trix elements for these heavy-light systems. 

One of the most important hadronic observables in B- 
physics is the Isgur-Wise form factor. Heavy quark sym- 
metry predicts that the matrix element of a heavy quark 
current operator in heavy-light mesons is described by a 
universal form factor 

(B{p')\b-f''b\B{p)) = Mb (v" + v'f') F{v ■ v') (3.1) 

Although this form factor does still depend on the quan- 
tum numbers of the light degrees of freedom, it becomes 
independent of the heavy quark mass as Mq —f oo. 
Knowledge of this form factor is very useful for deter- 
mining the CKM matrix element Vbc from semi-leptonic 
decays B D*ev^ because it allows one to extrapolate 



decay data to the zero recoil point v 



1. Within 



the context of the transverse lattice, the form factor is 
most conveniently extracted from the -|- component of 
the current operator, for which we find 



F{vv')=F'^'^\vv')+F'^^\vv'), 



where 



F(2)(:r) 
F(3)(x) 



2~x 
2 
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dwil>s{z, w)'ip 



z — X w 



1 — X ' 1 — X 



(3.2) 

(3.3) 
(3.4) 
(3.5) 



are the contributions from the two and three particle 
Fock components respectively and refers to the sum- 
mation over the discrete spin and orientation label for 
the meson. The velocity transfer is related to light- 
cone plus momentum transfer x through the relation 



V ■ V 
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I . Note that the normalization 

l — x J 

condition F{v-v' = 1) = 1 is automatic in this framework 
because of the normalization condition for the wave func- 
tion. Note also that terms that are off-diagonal in Fock 
space, and which in general have to be included when 
calculating current matrix elements when = — p^ 



oo (at 



is nonzero fl^ , are suppressed in the limit Mb 
least when considering the 'good' component). 

In the numerical calculations we evaluated Eq. (|3.2| ) 
for a variety of heavy quark masses and extrapolated to 
Mb oo. Results for the extrapolated form factor are 
depicted in Fig. |l|, where we multiplied the form factor 
with Vbc = 0.0424 [^| in order to be able to compare 
to experimental data |1^. Overall, our calculated form 
factor is in fairly good agreement with the data, although 



the experimental error bars are quite large. We should 
perhaps point out that if we include the data at larger 
values oiv-v' into account then we would obtain a slightly 
better fit to the data for values of Vbc in the range Vbc — 
0.0415 - 0.0420. The functional form of the fit to our 
calculated form factor as presented in Fig. |l| is: 




FIG. 1. Isgur-Wise form factor (3^) in the heavy quark 
limit. Uncertainties from the Mb oo limit were very small 
for this observable. The experimental data is from Ref. [14]. 

F{v.v') = 1 - 1.398{v.v' - 1) + 1.722(i;.t;' - i f 

-0.203{v.v' - if - 1.824(i;.w' - 1)^ (3.6) 

Although no experimental data exists on 



(3.7) 



it would be very useful to have predictions for this ob- 
servable because of its role in B-meson phenomenology. 

Up to an overall scale factor which involves the _L lat- 
tice spacing, Jb can be conveniently expressed in terms 
of the two particle wave function on the transverse lattice 
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(3.8) 



We find (in the heavy quark limit) for a value of Mb 
5GeV and a± « 0.5/m O] 



/s = 238±20MeV 



(3.9) 



where the error bar reflects only the uncertainty in our 
numerical extrapolation to the heavy quark lim it. A 
much larger systematic uncertainty in Eq. (3^) arises 
due to our Fock space truncation. We expect that in- 
cluding more higher Fock components will significantly 
reduce the 2-particle Fock space amplitude and hence 
also the decay constant, but it is very difficult to esti- 
mate the size of the reduction. 
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In addition to the decay constant, we also studied the 
hght-cone distribution amphtude itself. In the LF ap- 
proach, the twist-2 distribution amplitude is given by 
the wave function in the two particle Fock space sector 
with the quark and the antiquark at the same transverse 
position. On the transverse lattice, the distribution am- 
plitude, which plays a role in large momentum transfer 
processes, is thus conveniently expressed in terms of the 
two particle Fock component wavefunction 



{x)^NiP{x), 



(3.10) 



where A'^ is a normalization constant. The scaling behav- 
ior of (b(x) is such that 



^{x) = lim 



1 



(3.11) 



scales in the heavy quark limit. Numerical results for 
4ioo{x) are shown in Fig. ^. Although a precise determi- 
nation of the power law behavior near x — > was 



the bare (unscreened) longitudinal string tension is very 
large. We expect that including higher Fock states will 
lead to a significant lowering of the longitudinal string 
tension and hence also A. 

The light-cone momentum distribution of the heavy 
quark inside a B-meson is relevant for the invariant mass 
distribution in inclusive decays of heavy-light systems be- 
cause it determines the amount of energy that is avail- 
able for the weak decay products of the heavy quark. Of 
course, to leading order in 1/Mf, the h quark carries 100% 
of the B-meson's momentum in the infinite momentum 
frame. We thus consider the first correction to this trivial 
result, i.e. 



/,(^)= hm Ml-T^ 

Mi,-+oo \ Alt, 



(3.14) 



where h{x) is the light-cone momentum distribution of 
the 6-quark. Results are shown in Fig. ||. Since the 2 
particle Fock component is the dominant component for 
i3-mesons in our calculations, fb{z) is to a good approx- 
imation described by \'4}ao{z)^ (Fig- §)• 
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FIG. 2. B-meson distribution amplitude in the heavy quark 
limit. Systematic uncertainties arising from the numerical 
extrapolation are indicated by the shaded area. 

not possible, we found that our numerical results are 
consistent with a power close to 0.5 below x—1 GeV. For 
the lowest moments, we find 



dz(j)ao{z) 



1.509 Gel^- 



Ml = ^° , !; ; , = 1.223 GeV. 

Jo azcpooiz) 



(3.12) 



(3.13) 



From Eqs^3.12D and (3.13), it is clear that A = |Mi 



0.9 GeV pr in our calculations turns out to be very 
large. This result is consistent with the 'binding ener- 
gies' Mb — Mb that we obtain and reflect the fact that 
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FIG. 3. Light-cone momentum distribution (3.14). The 
shaded area reflects uncertainties arising from the numerical 
extrapolation. 



IV. SUMMARY AND OUTLOOK 

We have performed the first non-perturbative study 
of mesons containing one heavy and one light quark in 
light front QCD based on the transverse lattice fomula- 
tion. The parameters that appear in the Hamiltonian for 
a heavy-light system are the same that appear also in 
the Hamiltonian for light-light mesons. Since we deter- 
mined the latter in a previous work pl using constraints 
from rotational invariance, we were thus able to make 
parameter free predictions for a variety of heavy meson 
observables. 
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Specifically, we studied the Isgur-Wise form factor, de- 
cay constant, light-cone distribution amplitude and light- 
cone momentum distribution function for B mesons in 
the limit Mj, oo. Our results for the form factor are 
consistent with experimental data and favor a CKM ma- 
trix element Vbc in the range 0.0415 - 0.0420. 

Our main sources of systematic errors are the Fock 
space truncation (in combination with a truncation of 
the transverse lattice Hamiltonian Q), the truncation 
of the Hilbert space within each Fock component and 
the extrapolation to the —^ oo limit. The uncertain- 
ties from Hilbert space truncations are well under control 
since we used basis function techniques and extrapolated 
the results in the dimension of the Hilbert space. The 
Mh oo are also well under control since we studied 
a number of different Mh values. The main source of 
systematic errors in our calculations arise from the Fock 
space truncation and it is difficult to estimate the size 
of the corrections from higher Fock components without 
actually including them. Although we expect little cor- 
rections from higher Fock components for the form factor 
at small velocity transfers (the main effect of higher Fock 
component would be to provide the quarks with an addi- 
tional intrinsic form factor which would only be probed 
at large momentum transfers), higher Fock components 
probably have a large effect on the decay constant which 
is proportional to the amplitude to find the meson in the 
two particle Fock component. We did not include higher 
Fock components in this work because, for consistency, 
this requires inclusion of a number of new operators in 
the _L lattice Hamiltonian — a project that would have 
been beyond the intended scope of this one. Nevertheless, 
it is clear that one possible future extension of this work 
should be to take higher Fock components into account. 

There are many other natural extensions of this work. 
For example, it would be very useful to repeat this study 
for Bs mesons. While the ub meson calculations that we 
performed here could make use of light-quark parameters 
determined in previous works [pTp|] , such a study involv- 
ing s quarks would first require determination of hopping 
parameters for s quarks using techniques similar to the 
used in Refs. |n|,|. 

Although we performed the numerical calculations for 
a variety of heavy quark masses, we did not allow the 
heavy quarks to propagate in the transverse directions 
because this would have forced us to introduce (and de- 
termine) transverse hopping parameters. The advantage 
of this approach was that we were able to study the 
Mh —^ oo limit without introducing any new parame- 
ters. The disadvantage is that we cannot use the results 
to study the 1/Mf, to the heavy quark observables, which 
would be extremely useful for the analysis of experimen- 
tal data, especially those involving charm quarks. 
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